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Abstract 

Blood feeding is an integral process required for physiological functions and propagation of the malaria vector Anopheles. 
During blood feeding, presence of the malaria parasite, Plasmodium in the blood induces several host effector molecules 
including microRNAs which play important roles in the development and maturation of the parasite within the mosquito. 
The present study was undertaken to elucidate the dynamic expression of miRNAs during gonotrophic cycle and parasite 
development in Anopheles Stephens!. Using next generation sequencing technology, we identified 126 miRNAs of which 17 
were novel miRNAs. The miRNAs were further validated by northern hybridization and cloning. Blood feeding and 
parasitized blood feeding in the mosquitoes revealed regulation of 13 and 16 miRNAs respectively. Expression profiling of 
these miRNAs revealed that significant miRNAs were down-regulated upon parasitized blood feeding with a repertoire of 
miRNAs showing stage specific up-regulation. Expression profiles of significantly modulated miRNAs were further validated 
by real time PCR. Target prediction of regulated miRNAs revealed overlapping targeting by different miRNAs. These targets 
included several metabolic pathways including metabolic, redox homeostasis and protein processing machinery 
components. Our analysis revealed tight regulation of specific miRNAs post blood feeding and parasite infection in An. 
Stephens'!. Such regulated expression suggests possible role of these miRNAs during gonotrophic cycle in mosquito. Another 
set of miRNAs were also significantly regulated at 42 h and 5 days post infection indicating parasite stage-specific role of 
host miRNAs. This study will result in better understanding of the role of miRNAs during gonotrophic cycle and parasite 
development in mosquito and can probably facilitate in devising novel malaria control strategies at vector level. 
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Introduction 

Vector transmitted diseases are widespread in tropics and sub- 
tropics. Of several vector transmitted diseases, malaria accounts 
for more than 200 million cases with India contributing to almost 
70% of malarial cases in South East Asia region [1]. Malaria was 
brought on the verge of eradication by the success of malaria 
eradication program by World Health Organisation in 1950s. But 
1990s saw resurgence in malarial incidence due to insecticide 
resistance in vectors, drug resistance in parasites, changes in vector 
behaviour and lack of infrastructure to fight the disease. India has 
an intricate vector dynamics with six primary vectors of the 58 
Anopheles species reported in the country [2] . Of these, Anopheles 
stephensi is considered as an important vector with its distribution 
throughout the Middle East and South Asia region. 

Plasmodium passes through sexual and asexual stages of its 
development in mosquito and mammals respectively in order to 
complete its life cycle. Gametocytes transmitted to mammalian 
host by an infected mosquito develop into a motile ookinete which 
traverses through midgut epithelium and develop into an oocyst. 
Sporozoites burst out of oocyst, travel through hemolymph and 
invade salivary gland acinar cells. These sporozoites are then 



transmitted to a mammalian host upon next blood meal of 
mosquito. The several motile stages of parasite initiate an array of 
humoral and cellular immune responses in the mosquito. 
Activation of these pathways culminates in production of 
antimicrobial effector molecules, which function as either agonist 
or antagonist during maturation of Plasmodium parasite [3]. 
Previous studies have shown that blood feeding and Plasmodium 
infection leads to differential expression of genes in the mosquito 
[4,5]. Among the various factors known to regulate gene 
expression, microRNAs (miRNAs) have emerged as the most 
important regulatory molecules. MiRNAs are 22 nucleotides (nt) 
non-protein coding small RNAs shown to play role in multiple 
biological processes including defense response of host against 
invading pathogen [6] . Primary miRNAs are transcribed by RNA 
polymerase II [7] and is processed off its flanking regions in 
nucleus by an RNAse III enzyme, Drosha [8] giving rise to 
precursor miRNAs (pre-miRNA). Pre-miRNA is then exported to 
the cytoplasm by exportin -5 RanGTP complex [9-10] where it is 
processed by Dicer, into 22 nt miRNA duplex [11-13]. One 
strand of miRNA duplex binds to 3 ' untranslated region (UTR) of 
target genes and regulates their expression post-transcriptionally. 
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A single miRNA has been shown to regulate expression of 
thousands of mRNAs involved in diverse biological pathways [14—15]. 

Recent efforts towards control of spread of malaria have 
focussed on blocking Plasmodium cycle within the mosquito [16— 
1 7] . An increased understanding of vector-parasite interaction is 
essential to identify vector specific molecules that could be 
manipulated to block maturation of Plasmodium. Several studies 
have been conducted to understand role of miRNAs during 
growth and development in Drosophila [18-20]. Small RNA 
sequencing and computational methods have been successfully 
applied to identify putative miRNAs and their targets in Drosophila 
[21-24] and non-Drosophiloid insects. In comparison, scant 
information is available on differentially expressed miRNAs in 
mosquito at different stages of parasite maturation [6] . A recent 
study on An. stephensi has identified 27 miRNAs by direct cloning 
method. Further, they identified miRNAs implicated in mosquito 
reproduction and longevity [25] . 

Our present study was conducted to understand the role of 
Anopheles miRNAs during Plasmodium infection. For this purpose, 
we performed high throughput sequencing on An. stephensi 
mosquitoes after blood feed and at different time points of 
Plasmodium development in the vector. We identified 126 miRNAs 
of which 17 were novel miRNAs. Further analysis of these 
miRNAs revealed regulation of 1 3 miRNAs during blood feeding 
and 16 miRNAs upon Plasmodium infection. Additionally, we 
performed target prediction of these miRNA clusters to deduce 
their potential role. Our bio-informatic analysis predicted these 
regulated miRNAs to be playing roles in pathways directly 
affecting reproduction and parasite development in mosquito. 

Materials and Methods 

Ethics statement 

Animal experiments were performed in accordance with 
National animal ethics guidelines of the Government of India 
after approval by Institutional Animal Ethics Committees of 
International Centre for Genetic Engineering & Biotechnology, 
New Delhi (Permit number: ICGEB/AH/201 1/01/IR-8). 

Mosquito rearing and infection 

An. stephensi from laboratory stock were reared under controlled 
conditions at 28+2°C, 70-75% humidity. Adult mosquitoes were 
fed with 2% sterile glucose solution and water soaked raisin. Stock 
solution of Plasmodium vinckei petteri 279 BY was thawed and injected 
into 5 weeks old BALB/C mice. When gametocytemia reached 
0.05%, naive 4—5 days old female mosquitoes were fed on infected 
mice. Midguts of mosquitoes on 5 th day post infection (dpi) were 
dissected, stained in mercurochrome and observed for the 
presence of oocysts under microscope. Presence of oocysts 
confirmed infection of mosquitoes by Plasmodium parasite. 

Sample preparation 

Samples were prepared from An. stephensi at five different 
conditions. Sugar fed adult female mosquitoes 5-6 days old (SF) 
were collected as control sample. Blood fed (BF 42 h and BF 5d) 
and parasite infected mosquitoes (iBF 42 h and iBF 5d) were 
collected at two time points, at 42 hours and five days post blood 
feeding. Whole body of all mosquito samples were stored in Trizol 
(Invitrogen) until RNA extraction. Total RNA enriched in small 
RNA population was extracted using miRNeasy kit (Qiagen) as 
per the manual's protocol. Quality and quantity of RNA was 
checked by using Agilent 2100 Bioanalyzer RNA Nano 6000 kit. 



Small RNA sequencing 

Illumina Truseq small RNA libraries were made as per the 
manufacturer's instructions (Illumina Inc). 1 \ig of total RNA was 
ligated with 3' and 5' adaptors followed by reverse transcription 
using RT primers. Following PCR amplification of the adaptor 
enriched fragments, amplified products were ran and small RNA 
population within length 140-160 base pairs (bps) were eluted 
from 6% TBE PAGE gel. Eluted product was precipitated using 
sodium acetate and ethanol and was dissolved in RNAase free 
water. These small RNA libraries were then sequenced using 
Illumina Genome Analyzer II. 

Computational analysis of small RNA sequencing data 

Raw data generated post Deep sequencing was processed 
following an in-house pipeline described in detail [26] with slight 
modifications. Briefly, mature and pre-miRNA sequences of 
available eight insect species namely Anopheles gambiae, Aedes aegypti, 
Culex, Drosophila melanogaster, Bombyx mori, Apis mellifera, Acyrthosiphon 
pisum, Tribolium castaneum were downloaded from mirbase database 
v. 19 [27]. To reduce data redundancy, 100% similar sequences 
were pooled together using CDHIT web tool [28]. All available 
non-coding RNA sequences (ncRNAs) and protein coding region 
of Ae. aegypti, An. stephensi and An. gambiae were downloaded using 
ncRNA database [29] and vectorbase [30] respectively. All 
downloaded sequences were indexed to create four different 
databases using Bowtie [31]. Data generated after sequencing all 
five libraries were analysed separately. 

Reads derived from deep sequencing were trimmed and filtered 
to fetch sequences having length > = 18 bases (nt). In each library, 
same sequences were pooled together to generate expression files 
(fasta) with their read count and unique ID used for further 
analysis. These unique reads were aligned against mature miRNA 
sequence database using Bowtie with zero mismatches and taking 
other parameters as default. 

For identification of novel miRNAs, unmatched sequences after 
known miRNA prediction were mapped to pre-miRNA sequenc- 
es. Further, to remove non-coding RNAs, sequences left 
unmatched to pre-miRNAs were mapped to ncRNA database. 
Sequences that remained unmatched after this step were mapped 
to coding region of Ae. aegypti, An. stephensi and An. gambiae to filter 
out the sequences falling in coding regions. The final sets of 
unmatched sequences were then matched to An.stephensi genome. 
These matched sequences were subjected to novel miRNA 
prediction. 75 nt flanking region from both sides of genome 
sequence that matched with small RNA sequence was fetched as 
precursor sequence. These precursor sequences were folded to 
generate secondary structures and their folding energies were 
calculated using RNAfold [32] and RNAplot [32]. Precursor 
sequences forming hairpin loop structure with energy < = 
— 20 KJ/mol and small RNA sequence lying in one arm of 
precursor were reported as novel Pre-miRNAs and mapping small 
RNA sequence as mature miRNAs. 

To analyze relative abundance and expression profiling of 
miRNAs, tags per million of total RNA reads (TPM) for each 
miRNA in all five libraries were calculated. TPM was compared 
between the libraries to identify miRNAs differentially expressed 
post blood feeding and infection in mosquito. 

Expression profiling of miRNA by real time PCR 

The Custom miRNA locked nucleic acid PCR primer sets 
(Exiqon) were designed using target miRNA sequences from 
closely related species. 10 ng of Total RNA was reverse 
transcribed into cDNA using universal cDNA synthesis kit 
(Exiqon). Real time PCR was set up with two biological replicates 
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each in triplicates using SYBR green master mix (Exiqon) 
following manufacturer's instructions in ABI one step detection 
system. 5.8 s rRNA was used as an endogenous control for 
miRNA expression profiling. Expression levels were then calcu- 
lated against SF as a calibrator using 2 AAC X method. 



binding energy < = — 20 Kcal/mol were selected. KOBAS 
analysis was performed to identify significant pathways targeted 
by individual miRNAs keeping the threshold pvalue^O.05 [34], 
These selected targets were further used to generate miRNA: 
mRNA interaction network, visualized by Cytoscape [35] . 



Northern blot 

Northern hybridizations were conducted using digoxigenin- 
labelled antisense miRCURY LNA probes (Exiqon). Total RNA 
enriched in small RNA population was extracted using mirVana 
miRNA isolation kit (Invitrogen) from female mosquitoes collected 
at 42 h post blood feeding. 1 0 |Xg of total RNA was loaded in a 
15% denaturing polyacrylamide gels. Gels were stained with 
ethidium bromide for verification of RNA quality before 
transferring to a nylon membrane. Following l-ethyl-3-(3-dimethy- 
laminopropyl) carbodiimide crosslinking, membranes were pre- 
hybridized in a rotating hybridization oven for 30 min at 37°C in 
hybridization buffer. 0.5 nM miRCURY LNA probe (Exiqon) was 
added to the same buffer at 37°C for overnight hybridization. All 
probes were designed against Drosophila miRNA sequences. The 
membranes were washed twice for 10 min each in a low stringent 
buffer (2X saline sodium citrate, 0. 1% sodium dodecyl sulphate) 
at RT and then once in washing buffer (IX SSC) at RT. The 
membranes were incubated for 3 h in blocking buffer (Roche) 
followed by 30 min incubation with Anti-DIG-alkaline phospha- 
tase fAb (Roche) in blocking buffer. The membranes were then 
washed in DIG washing buffer at RT for 15 min each and then 
incubated for 5 min in development buffer. CSPD substrate (1:100 
diluted in development buffer) was applied on to the membranes 
and incubated in dark for 10 min. Chemiluminescence signal was 
then measured in a Fluorchem machine (Protein Simple) to detect 
miRNA on the membrane. 

Statistical Analysis 

Statistical tests for identifying significant differentially expressed 
miRNAs were performed using edgeR module with few modifi- 
cations in the script. The p value cutoff was performed on the data 
with the significance threshold selected as 0.05. 

miRNA target prediction and pathway analysis 

Target prediction and pathway analysis were carried out using 
protocol reported elsewhere [26]. Briefly, to understand miRNA 
function in mosquito, we predicted their messenger RNA (mRNA) 
targets using RNA hybrid tool [33]. 3'UTR of An. gambiae 
downloaded from vectorbase were used for the prediction. The 
targets showing complementarity with miRNA seed region and 



Results 

Small RNA sequencing using Illumina Truseq chemistry was 
performed on An. stephensi mosquito samples under five different 
conditions. Small RNA libraries enriched for miRNAs were 
constructed following manufacturer's instructions from sugar fed 
naive female mosquitoes 5-6 days old (SF), female mosquitoes at 
42 hours (BF 42 h) and 5 days post blood feeding (BF 5d) and from 
female mosquitoes at 42 hours (iBF 42 h) and 5 days post infected 
blood feeding (iBF 5d). Post library quantification and cluster 
generation, the libraries were sequenced using sequencing by 
synthesis technology, (Illumina Inc.) and the raw reads analysed 
further using an in-house pipeline developed in our laboratory 
[26]. 

Complexity of small non coding RNAs between libraries 
as revealed by deep sequencing 

Sequencing of all five libraries yielded 2.5 xlO 7 reads. After 
trimming the adaptors, 2.38 x 10 7 (94%) of total reads were utilized 
for further processing (Table 1). A total of 4.4 xlO 6 , 4.5 xlO 6 , 
3.9x10 s , 4.7 x 10 s , 6.2 xlO 6 sequences from SF, BF 42 h, BF 5d, 
iBF 42 h and iBF 5d libraries respectively were utilized for in silico 
analysis (Figure 1). Among the libraries size distribution of reads 
varied between 12-35 nt. We found bimodal distribution with 
peaks at 20-23 nt and again at 32-35 nt in all libraries. 20-23 nt 
peak corresponds to miRNAs, while 32-35 nt peak represents 
longer piRNA-like small RNAs. 32-35 nt small RNAs are twice as 
abundant in mosquitoes at 5 days post infected blood feeding 
suggesting a possible role of longer sized small RNAs in mosquito 
during parasite maturation. Whereas, differently sized small RNAs 
are equally distributed in blood fed and infected blood fed 
mosquitoes at 42 h post feeding (Figure 1). 

Small RNA reads between 18-30 nt, accounting for 78% of 
total reads were selected for further analysis. After removing 
redundant sequences, a total of 2.4xl0 5 , 9.5xl0 5 , 12xl0 5 , 
5.5x10 , 17xl0 5 unique sequences were found in SF, BF 42 h, 
iBF 42 h, BF 5d and iBF 5d libraries respectively. Total of 
2.3 x 10 6 reads gave hit with database containing ncRNAs (rRNA, 
tRNA, piRNA, snoRNA) from all organisms. Remaining sequenc- 
es were then analysed for known and novel miRNA candidates 
(Table 1). 



Table 1. Composition of small RNA in sugar fed naive female mosquitoes 6-8 days old (SF), blood fed female mosquitoes at 
42 hours (BF 42 h), 5 days post blood feeding (BF 5d), female mosquitoes at 42 hours (iBF 42 h) and 5 days post infected blood 
feeding (iBF 5d). 







Reads post 






Reads mapped to 


Unique small 


Sample Name 


Raw reads 


adaptor trimming 


18-30 nt reads 


ncRNA reads 


known miRNA (%) 


RNA reads 


SF 


5478556 


4426195 


2883395 


199021 


786128 (21) 


244917 


BF 42 h 


4635098 


4554220 


4132388 


410203 


769570 (18) 


954582 


iBF 42 h 


4862848 


4795865 


4150299 


817275 


257887 (6) 


1267813 


BF 5d 


4000410 


3908552 


3489024 


223131 


1020270 (29) 


558168 


iBF 5d 


6434830 


6207543 


5290072 


651297 


283816 (5) 


1 742904 


TOTAL 


2541 1 742 


23892374 


19945178 


2300927 


3117671 


4768384 



doi:1 0.1 371 /journal.pone.0098402.t001 
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Figure 1 . Length distribution of small RNA reads in (A) sugar fed naive female mosquitoes 6-8 days old (SF). (B) Female mosquitoes at 
42 h post blood feeding (BF 42 h). (C) 5 days post blood feeding (BF 5d). (D) Female mosquitoes at 42 h (iBF 42 h) and (E) 5 days post infected blood 
feeding (iBF 5d). X axis represents small RNA read lengths in base pairs while Y axis represents number of reads. (F) Distribution of abundant and rare 
miRNAs in all samples. Blue and Red bars represents abundant and rare miRNAs respectively. 
doi:1 0.1 371 /journal.pone.0098402.g001 



Identification of known and novel miRNAs in An. 
Stephens'! by Deep sequencing 

We identified a total of 109 known miRNAs in all five libraries 
with 1 00 % match with known miRNAs from eight insect species of 
which 103 miRNAs were conserved in mosquitoes (Table 2). A 
relative comparative distribution of these 103 mosquito conserved 
miRNAs is as follows: 65, An. gambiae; 87, Ae. aegypti and 67, Culex. 
Twelve miRNAs were identified in all eight insect species used in 
our analysis, illustrating evolutionary conservation of miRNAs 
across Arthropods (Table 2). Six miRNAs (miR-2796-5p, miR- 
2796-3p, miR-2779, miR-133-5p, miR-190-3p, miR-iab-8) iden- 
tified in our study had homologues with mature miRNAs from 
insect species other than mosquitoes (Table 2). 

High throughput small RNA sequencing data analysis identified 
several novel miRNAs in An. stephensi. Using our pipeline for novel 
miRNAs identification, 1 7 small RNA sequences showed all 
characteristic signatures of a miRNA (Figure 2 and Figure SI) as 
described earlier [36] . These miRNA sequences have been named 
in a specific order as detailed in Table 3. miRNAs found on 5' and 
3' arms of the precursors were annotated as -5p and -3p 
respectively. These sequences were further tested for hairpin 



formation and seven sequences showed complementary 5p and 3p 
miRNAs elevating these sequences as bonafide miRNAs (Table 3). 
NCBI BLAST was performed with these sequences and the results 
showed high conservation of these sequences across species. In 
order to categorize the novel miRNAs into known miRNA 
families, the sequences were blasted against mirBase, a database 
containing sequences of known miRNAs. Six novel miRNAs were 
tentatively classified into known miRNAs families as they showed 
highly conserved seed region with known miRNAs from other 
species. Sequence extending beyond seed region also showed high 
level of conservation with known miRNAs (Figure 3). 

Identification of miRNAs with sequence mismatches 

In the present study, identification of An. stephensi miRNAs was 
carried out using small RNA reads with 100% identity with known 
mature miRNAs. In addition, we performed the analysis with up 
to two mismatches in the mature miRNA sequences. By this 
method, we identified the following eight miRNAs (ast-miR-2941- 
3p, ast-miR-263a-3p, ast-miR-312-3p, ast-miR-316-5p, ast-miR- 
317-5p, ast-miR-956-3p, ast-miR-971-3p, and ast-miR-3843-3p) 
(Table SI). 
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A previous study has identified An. stephensi miRNAs using An. 
gambiae genome dataset [25]. In the present investigation, we 
mapped mature miRNAs of An. .stephensi identified in this study 
against the genome of An. stephensi, to validate their presence in the 
mosquito. All 109 miRNAs reported here mapped perfectly with 
the genome. However, keeping in mind the low coverage of An. 
stephensi genome and possible inaccuracy in genome sequence, we 
. reduced the stringency of our mapping and allowed up to two 

g mismatches and conducted the analysis. We identified three 

| miRNAs (ast-miR-2940-3p, ast-miR-2944b-5p, ast-miR-2951-3p) 

5 containing up to 2 mismatches with genome of An. stephensi (Table 
I SI), although they show perfect match with mature miRNAs in the 
:§ current miRBase database [27]. 

ii 

3 Validation of known and novel miRNAs in An. stephensi 

s by Northern blot 

■■^ Having identified miRNAs from the five libraries, efforts were 

g taken to validate some of these miRNAs using Northern 

.9- hybridisation. Total RNA from BF 42 h was probed using locked 

£ nucleic acid (LNA) probes (Exiqon Inc) for miR-286, miR-309, 

6 miR-989, selected at random (Figure 4a). 

ii Similarly, candidate sequences showing miRNA characteristics 

m were annotated as novel miRNAs and further validated by 
Northern hybridisation. From our dataset of predicted miRNAs, 

^ five miRNAs were selected in random and Northern hybridisation 

^ was performed. Two miRNAs (ast-nv-3-5p, ast-nv-4-5p) possessed 

<c TPM (Tags per million) >500 and were detected by Northern 

^ hybridization (Figure 4b). The other three miRNAs (ast-nv-2-5p, 

| ast-nv-6-3p ast-nv-7-3p) were not detected probably due to low 

g expression in our experimental samples. 
£ 

I miRNA abundance in sugar fed, blood fed and infected 
« An. stephensi 

Distribution and relative abundance of different Anopheles 

J miRNAs under three physiological conditions ie, sugar fed, blood 

jg fed and parasite infected blood fed was computed. Total number 

S of known miRNAs varied between the libraries of SF, BF and iBF 

g samples. Maximum number of known miRNAs (n=107) were 

| found in BF 5d while BF 42 h and iBF 5d possessed the least 

ig number of miRNAs each (n = 100) (Table 2). Total of 102 and 103 

§■ known miRNAs were identified in SF and iBF 42 h respectively 

| (Table 2). Out of 109 known miRNAs, 90 were found in all five 

ii libraries (Table 2). miRNAs were classified as "abundant" or 

| "rare" based on the tags per million (TPM) value. TPM >1000 

Tl were classified as abundant while those with TPM< 1 0 were 

■3 classified as rare. In case of SF, 20 miRNAs were classified as 

II abundant while 27 were classified under rare category (Figure IF), 
g. One miRNA, miR-2796-5p was seen exclusively in SF (Table 2). 
.!_- ^ However, the number of reads for this miRNA was very less (<5 
§; S reads) and requires further validation. At BF 42 h, 18 and 28 
0 8 3 miRNAs were found to be abundant and rare respectively. Three 
■g g. £ g (miR-286b, miR-2944a-3p, miR-309) were expressed significandy 

tJ ~ § in BF 42 h with no reads present in SF, indicating blood meal 

» c|? induced expression of these miRNAs (Table 2). Similarly, 26 were 

™ £ S o abundant and 31 were rare in BF 5d (Figure IF). miR-315-3p was 

J i; '>> c exclusively expressed in BF 5d with read count less than five 

III! ( Table2 )- 

S j J E Furthermore, abundance of individual miRNAs also varied 

|j ^ between sugar fed, blood fed and infected female mosquitoes as 

c = Pj reflected by the percentage of miRNA reads to the small RNA 

11 z z 5 population in the libraries (Table 1). Our analysis showed 

Si £ | '5 significant down-regulation of many miRNAs upon Plasmodium 
infection (iBF 5d). Four miRNAs, namely, bantam- 3p, miR-8-3p, 
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ast-nv-l-3p 



ast-nv-2-5p 
ast-nv-2-3p 
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ast-nv-5-3p 
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Figure 2. Predicted secondary structures of Novel miRNAs. Precursor sequences were folded using RNA fold. Coloured sequences (Pink = 

— 5p; Red = — 3p) represent mature miRNAs on both arms of the precursor miRNA. 

doi:10.1371/journal.pone.0098402.g002 



miR-184, miR-281-5p were the most abundant in all five libraries 
(Table 2). miR-263a-5p showed maximum number of reads in BF 
42 h and iBF 42 h, miR-bantam-3p was most the abundant 
miRNA in SF, miR-281-5p in BF 5d and iBF 5d. 

We found relatively less number of miRNAs in Plasmodium 
infected mosquitoes, 13 at iBF 42 h and nine in iBF 5d that were 
abundantly expressed indicating overall down-regulation of 
miRNA expression post infection in mosquitoes. iBF 5d contained 
maximum number of miRNAs showing TPM<10 (n = 51) 
whereas iBF 42 h contained 38 such miRNAs (Figure IF). 

Among six miRNAs (miR-2796-5p, miR-2796-3p, miR-2779, 
miR-133-5p, miR-190-3p, miR-iab-8) conserved in the eight 
arthropod species but not conserved in mosquitoes, miR-2796-3p 
was moderately expressed in all samples with iBF 5d showing 
maximum expression. miR-190-3p was also moderately expressed 
in SF with no reads in BF42 h. Remaining four miRNAs showed 
very less abundance in all the stages (Table 2). 

Temporal regulation of miRNAs post blood feeding and 
infection in An. Stephens'! 

Exploiting the large data and deep coverage from high 
throughput sequencing of the libraries, we generated miRNAs 
expression profiles in blood fed and infected female mosquitoes. 



MiRNAs were profiled at two time points i.e. at 42 hours and 5 
days post infected blood feeding coinciding with late phase of 
midgut invasion by the parasite and initiation of sporozoites 
release from the oocysts. Using in-house developed Perl script, we 
calculated TPM for individual miRNAs and compared their 
expression profiles across different experimental conditions. 
EdgeR module was used to identify significantly regulated (p 
value^0.05) miRNAs between different experimental conditions. 
The final set of significantly regulated miRNAs was selected on the 
basis of their p value^0.05, fold change Si. 5 and TPM>10. 

With respect to blood feeding, 11 miRNAs (miR-2944a-5p, 
miR-92b, miR-989, miR-275, miR-281-3p, miR-281-5p, miR- 
306, miR-263a-5p, miR-7, miR-309 and miR-305-3p) were 
significandy up-regulated in BF 42 h compared with SF (Table 2, 
Figure 5). Whereas, we found only two miRNAs, namely, miR- 
190-3p and miR-124 down-regulated in BF 42 h compared with 
SF. In case of BF 5d, most of the miRNAs showed a similar profile 
with that of SF indicating the similarity in the physiological state of 
mosquito at this stage (Figure 5). We identified 15 miRNAs that 
were regulated in BF 5d compared with SF, possibly due to one 
round of blood digestion and egg production in BF 5d (Table 2). 

We noticed distinct variation in expression pattern of miRNAs 
in the following corresponding conditions, namely 42 hrs and 5 
days post Plasmodium infection. In order to understand Plasmodium 
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Table 3. List of novel miRNAs identified in An. Stephens:. 





S. No 


Length 


Sequence 


miRNA 


Tags Per Million (TPM) 












miRNAs on both arms of Pre-miRNA 




SF 


BF 42 h 


iBF 42 h 


BF 5d 


iBF 5d 


1 


22 


ATTAGAATGTGGAATCTGTTTT 


as-nv-1 -5p 


5.2 


2.5 


2.8 


7.9 


0.6 


2 


22 


AAAAGTTTTCATATTCTTGCGG 


as-nv-1 -3p 


3.2 


1 


0.8 


1.4 


0.3 


3 


22 


AAACAAGTCAAGGAATCAGGGA 


as-nv-2-5p 


30.8 


17.9 


6.9 


19.2 


5.9 


4 


21 


TTGATTTCTCGACTTGATTGT 


as-nv-2-3p 


16.2 


15.5 


5.5 


28.2 


2.1 


5 


28 


TGTAATGTACTCTCGTTCTGGACGATTT 


as-nv-3-5p 


4.7 


643 ¥ 


200.9 


1350.1 


177.3 


6 


21 


TTTCGGATATGAATCAAAGTA 


as-nv-3-3p 


5.2 


78.7 


35.7 


34.4 


12.1 


7 


22 


GTAAGAGAAGCAATCGCGGGTT 


as-nv-4-5p 


520 ¥ 


356.6 


133 


817.1 


93.3 


8 


23 


CGCGTGACTTGCTTCACTTTCGC 


as-nv-4-3p 


1 


1 


0.6 


2.7 


0.6 


9 


26 


AATGGTCTACTGCCACGTGCATCGCA 


as-nv-5-5p 


1.8 


11.6 


6.7 


19.7 


2.1 


10 


26 


TATGCGGAAACTTTTAAAAGGATGGG 


as-nv-5-3p 


9.3 


9.4 


2 


8.4 


1 


11 


25 


TGCACCCCCGGGCTGAGAAGAATCC 


as-nv-6-5p 


8 


17.6 


4.9 


39.7 


4.6 


12 


26 


TAAGATGGACAGCCGGGAAACTGATC 


as-nv-6-3p 


2.7 


24.1 


2.4 


19.2 


5.5 


13 


28 


TGCAAGAAGGAACTATACTCCGACGCCT 


as-nv-7-5p 


2.7 


2.8 


1 


4.2 


0.4 


14 


27 


TGTATTGGTCGTCAGCAATGTAGTCCT 


as-nv-7-3p 


56 


34.7 


9.2 


34.4 


3.5 


miRNAs on one arm of Pre-miRNA 


15 


23 


ACTTTGGTGATTTAGCTGTATGT 


as-nv-8 


181.7 


529.6 


155.8 


571.6 


193 


16 


23 


TGGTTTATCTAGTCTGTCGAAAC 


as-nv-9 


26.1 


14.8 


2.4 


26.2 


4.6 


17 


25 


TGATGTAATGTACTCTCGTTCTGGA 


as-nv-1 0 


10 


16.1 


8 


42.4 


6.3 


18 


27 


TATGACCAGAGGATAGATGTGACTACT 


as-nv-1 1 


8 


22.4 


120.7 


11.4 


141.8 


19 


23 


CGGGGGATGCTGTACGCTCTATG 


as-nv-1 2 


6.7 


16.1 


4.3 


10.7 


2.4 


20 


25 


TACGAGAAGATGGTCGCAAAAGCTC 


as-nv-1 3 


6.3 


14 


6.1 


15.9 


6.5 


21 


24 


ACGACGTATGAGGATACCCTGAAA 


as-nv-1 4 


5.4 


85.4 


7.1 


17.9 


6.3 


22 


25 


TCCTTCTGTTGACCTGGCGCTCGAC 


as-nv-1 5 


5.1 


26.7 


5.3 


21.7 


4.8 


23 


19 


TCACAGTCTTACTCGGGTG 


as-nv-1 6 


8 


8.1 


2.2 


13.4 


1.2 


24 


26 


TAGAAAGGCGCCGTAGAATCTTGGCC 


as-nv-1 7 


4.3 


8.1 


3.9 


12.7 


0.9 



¥ novel miRNAs validated by Northern hybridization. miRNA ast-nv-4-5p was validated using sample SF; miRNA ast-nv-3-5p was validated using sample BF 42 h. 
doi:1 0.1 371 /journal.pone.0098402.t003 



stage-specific regulation, the miRNAs were profiled at two specific 
conditions coinciding with effect of parasite first invasion in 
Anopheles midgut at 42 hrs and late phase of Plasmodium develop- 
ment in midgut [37]. It is interesting to note that distinct subsets of 
miRNAs showed regulation in the sets analysed namely iBF 42 h 
vs BF 42 h and iBF 5d vs iBF 42 h. In iBF 42 h, four miRNAs 
(miR-124, miR-137, miR-1000, miR-932) were significantly up- 
regulated compared to BF 42 h. However, none of the identified 
miRNAs were down-regulated in iBF 42 h compared to BF 42 h. 
In iBF 5d, four miRNAs (miR-1 1 75-3p, miR-1174, miR-281-3p 
and miR-281-5p) were significantly up-regulated, whereas 10 
miRNAs (miR-285, miR-2944a-5p, miR-309, miR-210-3p, miR- 
1891, miR-981, miR-315-5p, miR-932, miR-124 and miR-7) were 
significantly down-regulated when compared with iBF 42 h 
showing parasite stage specific expression in mosquito (Table 2, 
Figure 5). 

Expression profiling of novel miRNAs in blood fed and 
infected mosquito 

TPM values of novel miRNAs were compared between libraries 
to identify miRNAs that showed regulated expression post blood 
feeding and infection in mosquito. miRNAs ast-nv-3-3p, ast-nv-3- 
5p, ast-nv-5-5p, ast-nv-6-3p, ast-nv-8, ast-nv-12 and ast-nv-14 
were found to be up-regulated in BF 42 h (Figure 5B, Table 3). 



As observed in the case of known miRNAs, most of the novel 
miRNAs were down-regulated in Plasmodium infected mosquito 
samples. The miRNA ast-nv-1 1 was up-regulated in iBF 42 h and 
iBF 5d. Whereas, ast-nv-3-5p, ast-nv-3-3p and ast-nv-5-5p were 
significantly up-regulated in iBF 42 h compared with SF 
(Figure 5B, Table 3). ast-nv-6-3p was found to be up-regulated 
in iBF 5d compared with iBF 42 h. Our analysis also revealed 
three miRNAs, ast-nv-2-5p, ast-nv-1 -3p, ast-nv-7-3p to be 
significantly down-regulated in all stages of blood fed and infected 
mosquito compared with SF (Figure 5B, Table 3). 

Validation of miRNA profile by real time PCR 

To verify small RNA sequencing results described above, the 
expression profile of six differentially expressed miRNAs were 
validated by quantitative real time polymerase chain reaction 
(qRT-PCR). All six miRNAs showed similar expression profiles as 
revealed by small RNA sequencing analysis. Both small RNA 
sequencing and RT-PCR showed that miR-309 and miR-989 
were up-regulated in BF 42 h compared with SF (Figure 6A and 
6B). Similarly, miR-309 did not show significant expression 
difference in iBF 42 h and iBF 5d when compared to SF 
(Figure 6A) whereas miR-989 was highly up-regulated in both time 
points of infected mosquito compared to SF (Figure 6B). Both the 
miR-210 (Figure 6D) and miR-285 (Figure 6E) were exclusively 
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199 181.799729709799 996 ACTTTGGTGATTTAGCTGTATGT 



3557 6.7536044169303 37 CG G G G G ATG CTGTACG CTCTATG 



UserSeq 
aga-:tiR-9a 



: | ;uuu;;uj^uuua;;uguaug 22 

I I I I I I I II I I I I I I I I I I 
2 cuuugguuaucuagcuguaug 22 



UserSeq 
dps-ir.iR-965-5p 



; I z z z z z i jj z : u :ua : z cu cuaua 21 

III I I I I I I I I I I I I I I 
2 gggucaaacuguacguucuaug 23 



1048 26.1017684221901 143 TGGTTTATCTAGTCTGTCGAAAC 



5236 4.38071637854938 24 TAGAAAGGCGCCGTAGAATCTTGGCC 



UserSeq 

aae-rr.iR-2940-5p 



1 .£^IuajJ:uagucugucga 2 0 

I I I I I I I I I I I I I I I I I I 
1 ugguuuaucuuaucugucga 2 0 



UserSeq 
crir.-ir.iR-60-5p 



1 u hgaaaqql sqccquaqaauc 20 

I I I I I I I I I I II I I 1 1 
4 uagaaaggugacauaaaauc 23 



3070 8.03131336067387 44 TCACAGTCTTACTCG G GTG 



2666 9.30902230441744 51 TGCGATGCACGTGGCAGTAGACCATT 



UserSeq 
asu-nnR-43d-3p 



1 ± r^l .-acueaa 

I I I I I I I I I I I III 
3 ucacagucuuauccgg 



!6 UserSeq 

: = gga-ru.R-1564-5p 



2 I r::aur:i acquqqcaquaqa 21 

II III II I I I I I I I I I 
1 qcaauqaacacggcaguaga 20 



Figure 3. Family-wise clustering of novel miRNAs. Alignment results of novel miRNA sequences with known mature miRNAs present in 
mirbase. Seed region of Novel miRNAs are enclosed in a box. 
doi:1 0.1 371 /journal.pone.0098402.g003 



down-regulated at iBF 5d when compared to other experimental 
conditions. miR-932 was significandy up-regulated at iBF 42 h 
compared with BF 42 h whereas it was down-regulated at iBF 5d 
compared to iBF 42 h (Figure 6C). At iBF 5d, miR-1174 was 
identified to be significantly up-regulated compared to iBF 42 h 
(Figure 6F). All RT-PCR amplified products were cloned and 
sequenced to validate amplification of correct miRNAs. BLAST 
results revealed specific miRNAs amplified by RT-PCR, thereby 
further validating small RNA sequencing. 

miRNA target prediction and pathway analysis 

To understand the function of miRNAs, it is obligatory to study 
their targets. As miRNAs work in clusters to regulate various 
biological processes, we aimed to predict targets for all miRNAs 
that were regulated at different experimental conditions. We 
predicted targets for all those miRNAs that were significantly 
modulated, namely, 1 2 miRNAs in BF 42 h, four miRNAs in iBF 
42 h and 14 miRNAs that were regulated in iBF 5d. miR-190-3p 
was eliminated from target analysis as its reads were not found in 
BF 42 h though it was moderately expressed in SF and classified as 
down-regulated miRNAs at BF 42 h. Maximum number of 
mRNAs were targeted by miR-989 (n= 1222), miR-1000 (n= 741) 
and miR-281-3p (n = 284) that were up-regulated at BF 42 h, iBF 
42 h and iBF 5d respectively (Table S2). miR-124 was the only 
down-regulated miRNAs at BF 42 h that targeted (n = 505) 
mRNAs in mosquito (Table S2). Among down-regulated miRNAs 
at iBF 5d, maximum and minimum number of mRNAs were 
targeted by miR-124 (n = 505) and miR-7 (n = 99) respectively 
(Table S2). The targets were analysed by KOBAS and those falling 
under significant pathways (p value<0.05) were studied. miRNAs 
regulated at BF 42 h targeted 23 different pathways with eight of 



them being shared between two different miRNAs (Figure 7A). 
Two pathways, carbon metabolism and protein processing in 
endoplasmic reticulum were commonly targeted by down-regu- 
lated (miR-124) and up-regulated (miR-989) miRNAs. Post 
infection, we identified nine and 28 different pathways targeted 
by miRNAs regulated at iBF 42 h and iBF 5d respectively 
(Figure 7B and 7C). Among these pathways, two were shared 
between two or more miRNAs at iBF 42 h and ten different 
pathways by iBF 5d (Figure 7B and 7C). At iBF 5d, five pathways 
were commonly targeted by up- and down-regulated miRNAs. 
Oxidative phosphorylation was targeted by miR-1174 and miR- 
309, Proteasome was common between miR-1174, miR-1175-3p 
and miR-2944a-5p, RNA transport commonly targeted by miR- 
1175-3p, miR-981 and miR-1891, protein processing in endo- 
plasmic reticulum targeted by miR-1175-3p, miR-285, miR-124, 
miR-210-3p and miR-1891 and N-Glycan biosynthesis was 
targeted by miR-1175-3p and miR-1891 (Figure 7C). Among all 
pathways, we identified protein processing in endoplasmic 
reticulum to be regulated by maximum number of miRNAs 
namely miR-124, miR-1000 and miR-137 in iBF 42 h and by 
miR-1175-3p, miR-285, miR-124, miR-210-3p and miR-1891 in 
iBF 5d (Figure 7B and 7C). 

miRNA:mRNA interaction was analyzed through network 
generation as shown in Figure 8. Several transcripts were found 
to be targeted by two or more miRNAs. At BF 42 h eight 
transcripts (AGAP002499, AGAP006966, AGAP007642, 
AGAP011960, AGAP02728, AGAP004867, AGAP004362 and 
AGAP0 11374) were targeted by two or more miRNAs up- 
regulated at BF 42 h (Figure 8A). Whereas, 1 1 transcripts 
(AGAP010866, AGAP000184, AGAP001312, AGAP002192, 
AGAP004742, AGAP011551, AGAP002564, AGAP010557, 
AGAP005404, AGAP004192 and AGAP008327) were found to 
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Figure 4. Northern hybridization of An. Stephens'! miRNAs. (A) Validated presence of (1) miR-286, (2) miR-309 and (3) miR-989 in An. stephensi 
via DIG based Northern hybridization. (B) Northern hybridization of novel miRNAs, (4) ast-nv-3-5p, (5) ast-nv-4-5p in An. stephensi with TPM>500. 
10 ug of Total RNA from BF 42 h mosquito was used for northern blots of (1), (2), (3) and (4) whereas (5) was validated using total RNA from SF. RNA: 
Total RNA run on PAGE; L: Ladder. 
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be commonly targeted by up- and down-regulated miRNAs at BF 
42 h (Figure 8A). At iBF 42 h, two or more up-regulated miRNAs 
targeted 11 different mRNAs (AGAP000184, AGAP001986, 
AGAP002192, AGAP004192, AGAP007975, AGAP006569, 
AGAP008918, AGAP004098, AGAP004609, AGAP001759 and 
AGAP004176) (Figure 8B). Target analysis at iBF 5d revealed 13 
mRNAs (AGAP007538, AGAP011504, AGAP005535, 
AGAP008837, AGAP004867, AGAP000862, AGAP004020, 
AGAP001986, AGAP007001, AGAP005861, AGAP005366, 
AGAPO 11888 and AGAP009570) common between up- and 
down-regulated miRNAs (Figure 8C). Furthermore, we identified 
nine mRNAs (AGAP002499, AGAP000832, AGAP000747, 
AGAP003560, AGAP002728, AGAP002192, AGAP011551, 
AGAP004362 and AGAP007309) and one mRNA (AGAP000308) 
were targeted by up- and down-regulated miRNAs at iBF 5d 
respectively (Figure 8C). 

Discussion 

RNA interference is the most important and primary defence 
employed by insects to protect themselves from pathogens. Several 
components of this important phenomenon are involved in various 
biological processes, miRNAs being the most important of such 



molecules. Several miRNAs have been identified in various 
mosquito species and some of their roles in insect reproduction 
and infection have been explained [6,25]. However, their role in 
Plasmodium development is poorly understood. This study was 
conducted to identify and elucidate role of miRNAs post blood 
feeding and Plasmodium infection in An.stephensi. 

The present study is the first report where small RNA reads 
were mapped on to An. stephensi contigs for identification of novel 
miRNAs which has led to the identification to a large number of 
known and novel miRNAs. Seven pairs (-5p and -3p) of miRNAs 
were identified that fulfilled all criteria required to be classified as a 
miRNA [36]. Few were further classified into known miRNA 
families based on seed homology. Additionally, some of novel 
miRNAs identified lacked any sequence homology with known 
miRNAs, thereby identifying new families of miRNAs that might 
play important role during Plasmodium development in vector. 

To elucidate the role of miRNAs listed in this study, information 
related to their expression profile post blood feeding and infection 
is necessary. Here, expression profile of miRNAs at two time 
points post blood feeding and parasite infection were investigated 
using small RNA sequencing and further validated by qRT-PCR. 
In our study, we found significant down-regulation of miRNAs 
post infection of mosquito by Plasmodium parasite. Due to 
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Figure 5. Heat Map of differentially expressed miRNAs present in small RNA libraries of An. Stephens/. (A) Expression profile of miRNAs 
in sugar fed naive female mosquitoes 6-8 days old (SF), 42 h post blood fed female mosquito (BF 42 h), 42 h post infected blood fed female 
mosquito (iBF 42 h), 5 days post blood fed female mosquito (BF 5d), 5 days post infected blood fed female mosquito (iBF 5d). Colour gradation from 
light green to dark red indicates relative increase in miRNA expression. (B) Heat Map of differentially expressed novel miRNAs present in all small RNA 
libraries of /An. Stephens'!. Expression profile of miRNAs in sugar fed naive female mosquitoes 6-8 days old (SF), 42 h post blood fed female mosquito 
(BF 42 h), 42 h post infected blood fed female mosquito (iBF 42 h), 5 days post blood fed female mosquito (BF 5d), 5 days post infected blood fed 
female mosquito (iBF 5d). Colour gradation from light green to dark red represents relative increase in miRNA expression. 
doi:1 0.1 371 /journal.pone.0098402.g005 



co-evolution of host-parasite crosstalk, hosts are known to resist 
parasite infection whereas parasite manipulate host to increase its 
chances of survival and transmission rates [38]. Also, Plasmodium 
infection is known to alter transcriptome and proteome profile in 
mosquito [4,39]. Therefore, miRNA regulation in our study could 
have been brought about by targeting the components involved in 
host miRNA biogenesis pathway either by molecules present in 
blood or the invading parasite. Keeping in view the process of 
establishment of parasite in insect, we hypothesize parasite derived 
factors might play a role in the regulation of miRNAs involved in 
defense response of host against invasion. 

Blood feeding in hematophagus insects, such as mosquito is 
essential for reproduction and is exploited by Plasmodium parasite 
to complete its life cycle. Identification of miRNAs playing role in 



insect reproduction is crucial to control mosquito population and 
hence parasite transmission from one host to another. We 
identified few miRNAs that were up-regulated in mosquito post 
blood feeding. Such induced expression indicates towards their 
involvement in mosquito reproduction. In the present investiga- 
tion miR-275 was up-regulated upon blood meal ingestion. 
Employing anti-miRNA in feeding experiments, miR-275 was 
shown to play a direct role in blood digestion and egg development 
in Ae. aegypti [40]. The observed up-regulation of miR-275 in 
present set of experiments also suggests its role in blood digestion 
and egg development and correlates with the observed physiolog- 
ical situation. Another study on Aedes reported induced expression 
of aae-miR-375 in presence of blood meal [41]. We did not 
identify miR-375 regulated at BF 42 h when compared with SF, 
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Figure 6. Expression profiling by Real time PCR of up-regulated and down-regulated miRNAs. (A) miR-309,(B) miR-989,(C) miR-932, (D) 
miR-210,(E) miR-285 and (F) miR-1 174 were profiled in sugar fed naive female mosquitoes 6-8 days old (SF), female mosquitoes at 42 hours (BF 42 h) 
and 5 days post blood feeding (BF 5d) and from female mosquitoes at 42 hours (iBF 42 h) and 5 days post infected blood feeding (iBF 5d). Y axis 
depicts fold change in miRNA expression in samples compared with sugar fed naive female mosquitoes 6-8 days old (SF), taken as 1. 
doi:10.1371/journal.pone.0098402.g006 



probably due to insect specific regulation of this miRNA. In 
addition to these miRNAs, miR-309 was found to be induced by 
blood feeding in mosquito. This miRNA has been implicated in 
facilitating normal embryo development in Drosophila and may be 
involved in similar biological process in An. stephensi [22,23]. Also, 
miR-309 was not identified in SF in our analysis that correlates 
with previous study where it was identified to be involved in 
turnover of maternal mRNA during maternal to zygotic transition 
in fly [42]. Another miRNA, miR-989 was shown to be up- 
regulated post blood feeding in mosquitoes and was abundandy 
expressed in ovaries of blood fed mosquitoes [6,25] . Our analysis 
also revealed up-regulation of this miRNA in small RNA 
sequencing and RT-PCR analysis of BF 42 h mosquito. 

In this study, miRNAs were profiled at two time points i.e. at 
42 hours and 5 days post infected blood feeding coinciding with 
late phase of midgut invasion by the parasite and initiation of 
sporozoites release from the oocysts. Previous study [6] has shown 
that miR-1174 and miR-1 175 were down-regulated in midgut 
tissue at 24-48 hrs post infected blood feeding. Our study showed 
up-regulation of the above mentioned miRNAs in whole body of 
mosquitoes at iBF 5d. The same study also identified miR-989 



up-regulated in midgut of infected female mosquito. However, in 
our study, miR-989 was not significantly regulated at iBF 42 h 
whole body compared with BF 42 h. The reason for such 
discrepancy could be the different tissues, namely midgut used for 
profiling of these miRNAs in above mentioned study whereas our 
study was performed on whole body tissues. Additionally, two 
miRNAs (miR-281-5p and miR-281-3p) were up-regulated at iBF 
5d compared with iBF 42 h. These miRNAs may have parasite 
stage specific function through Toll pathway as miR-281 has 
putative targets involved in Toll pathway in D. melanogaster [43]. 
Same study identified multiple genes involved in immune 
pathways of D. melanogaster as putative targets for miR-281, miR- 
124, miR-315, miR-309 and miR-210 that were found to be 
regulated in mosquitoes during parasite maturation in this study. 

Numbers of genes are altered at transcriptome as well as 
proteome level during the course of blood feeding and parasite 
development in the mosquito [4,39]. Few of these altered 
molecules might serve as agonist or antagonist against Plasmodium 
development. A number of factors are known to regulate 
expression of such genes including miRNAs. Therefore, in order 
to link miRNA regulation observed in this study with mRNAs that 
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Figure 7. KOBAS analysis of miRNA targets predicted by RNA hybrid. Pie charts represents significant pathways (P value<0.05) targeted by 
miRNAs regulated at (A) BF42 h, (B) iBF 42 h and (C) iBF 5d respectively. Pie chart area for each pathway represents percentage of regulated miRNAs 
targeting it. miRNAs are listed in the pie area of respective pathway it is targeting. Pathway marked star are targeted by maximum number of 
miRNAs. 
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might play role during egg development and parasite life cycle, we 
carried out preliminary target prediction of selected miRNAs. 
Each miRNA targeted multiple transcripts. Many of these 
transcripts were targeted by two or more miRNAs. In mosquito, 
blood meal is required for initiation of egg maturation involving 
numerous pathways. One such pathway is TOR signalling 
pathway, known to activate production of yolk protein precursors 
for egg development in blood fed mosquito [44]. This pathway 
was found to be targeted by miR-305-5p, highlighting its 
importance in events occurring post blood feeding in mosquito. 
Blood feeding in mosquito also induces major metabolic challenge 
including change in redox homeostasis due to rapid blood 
digestion and assimilation of nutrients present in the blood meal. 
In our study, miRNAs regulated post blood meal were identified to 
target mRNAs involved in amino acid metabolic pathways 
required for oogenesis and preventing accumulation of amino 
acid metabolites in the mosquito. Number of pathways involved in 
maintaining redox state such as pentose phosphate pathway, 
citrate cycle and oxidative phosphorylation were found to be 
targeted by miR-124, miR-2944a-5p, miR-305-5p and miR-309. 
miR-309 was also identified to target nitrogen metabolic pathway 
essential for protein metabolism and detoxification of ammonia in 
the mosquito. Among the common targets, few mRNAs 



(AGAP004742, AGAPO 11551 and AGAP002564) coding for 
pyruvate carboxylase, protein transport protein SEC 31 and 
aldolase respectively were previously found to get regulated post 
blood feeding and embryonic development in mosquito [45-46]. 

As in response to blood feed, successful infection by Plasmodium 
in mosquito also triggers multiple changes in metabolic, redox and 
immune related processes. Previous studies have shown that 
metabolic and redox state of mosquito can alter mosquito 
susceptibility to Plasmodium infection [47-48]. At iBF 42 h, up- 
regulated miRNAs namely miR-124 and miR-932 were identified 
to target pentose phosphate pathway and citrate cycle respectively. 
These targets were probably required for the controlled ROS 
generation in mosquito for optimal Plasmodium development and 
transmission to mammalian host. Further, we identified miR-137 
targeting mRNAs playing role in TGF-beta signalling pathway, 
essential for mosquito immunity against Plasmodium parasite. One 
miRNA, miR-124 is of utmost importance as it was significandy 
down-regulated at BF 42 h whereas was up-regulated at iBF 42 h. 
It was identified to target 25 mRNAs lying in carbon metabolism, 
pentose phosphate pathway and protein processing in endoplasmic 
reticulum pathways. At iBF 5d, miR-281-3p was found to target 
mRNAs in TGF-beta signalling pathway. Glutathione metabolic 
pathway is a major redox pathway, whose genes were found to be 
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Figure 8. miRNA-mRNA interaction network. Interaction network of regulated miRNAs and their targets at (A) BF 42 h, (B) iBF 42 h and (C) iBF 

5d. miRNAs are diamond shaped whereas targets are circular. Up-regulated miRNAs are highlighted in red whereas down-regulated are highlighted 
in green colour. Transcripts targeted by two or more miRNAs are marked in yellow colour. Significance of the Transcripts marked (*) are explained in 
the discussion. 
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targeted by miR-285. Further, many more metabolic pathways 
were identified to be targeted by miR-1174, miR-2944a-5p and 
miR-309 regulated at iBF 5d. Analysis of common targets revealed 
some important transcripts such as AGAP004192, AGAP000747 
and AGAP01 1888 coding for heat shock protein, MAPKKK and 
exportin-5 respectively. Previous studies have shown their 
importance in pathogenesis of infectious diseases. Hosts are 
known to increase heat shock protein production in response to 
pathogen encounter [49]. Similarly, MAP protein kinases are a 
part of MAP Kinase signalling pathway playing role in anti- 
Plasmodium defence in mosquitoes [50] . Exportin -5 is involved in 
miRNA biogenesis pathway in different organisms. Such compo- 
nents of miRNA biogenesis pathway are known to get regulated in 
host cells to enhance the survival of infectious agents [51-52]. In 
our study, exportin-5 was found to be targeted by miR-1175-3p 
and miR-1891. miR-1175 was shown to get up-regulated post 



blood feeding whereas infection with Plasmodium parasite causes its 
down-regulation in the midgut of female mosquito [6]. A previous 
study has shown role of miR- 1891 in longevity and fecundity of 
adult Aedes albopictus mosquitoes [53]. 

Parasites are known to alter host responses for their optimal growth 
and development. Infectious agents like HIV-1 regulate miRNA 
silencing pathway of the host during its replication [54]. Bacterial 
effector proteins are also known to suppress miRNA pathway to 
cause disease [55] . Uninfected blood is known to harbour number of 
factors resulting in differential expression of miRNAs. Differential 
expression of miRNAs upon parasite infection can either be a result 
of defence response of host against invading parasite or brought 
about by parasite itself to sustain its development. Therefore, whether 
the differential expression of miRNAs found in this study is attributed 
to host responses or are caused by invading parasite effector 
molecules needs further investigation. 
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This is the first comprehensive study to profile a large number of 
An. stephensi miRNAs at different time points post blood feeding 
and Plasmodium infection. Differentially expressed miRNA in BF 
and iBF were identified by small RNA sequencing. These were 
further validated by real time PCR. Novel miRNAs also showed 
regulated expression in BF and iBF. Many miRNAs identified 
were conserved in insects although few were identified for the first 
time in Anopheles. Target prediction revealed important mRNA 
candidates mapped to diverse biological pathways. miR- 
NA:mRNA interaction network point towards complex regulation 
of processes post blood feeding and infection in the vector. 
Elucidation of regulated miRNA functions is an important 
challenge that will provide strong foundation for understanding 
parasite development and mosquito-parasite interactions. 

Supporting Information 

Figure SI Figure SI is a figure showing predicted 
secondary structure of novel miRNAs. Mature miRNA 
were found only on one arm of pre-miRNA, represented by 
coloured sequences. 
(TIF) 

References 

1 . Dash AP, Valecha N, Anvikar AR, Kumar A (2008) Malaria in India: challenges 
and opportunities. Journal of Biosciences 33:583—592. 

2. Nagpal BN, Sharma VP (1995) Indian Anophelines. Science Pub-ushers 1-416. 
.3. Osta MA, Christophidcs GK, Kafatos FC (2004) Effects of Mosquito Genes on 

Plasmodium Development. Science 303:2030—2032. 

4. Felix RC, Miiller P, Ribeiro V, RansonH SilvciraH (2010) Plasmodium infection 
alters Anopheles geunbiae detoxification gene expression. BMC Genomics 1 1:312. 

5. Dana AN, Hillenmeyer ME, Lobo NF, Kern MK, Romans PA, ct al. (2006) 
Differential gene expression in abdomens of the malaria vector mosquito, 
Anopheles gambiac, after sugar feeding, blood feeding and plasmodium berghei 
infection. BMC Genomics 7:119. 

6. Winter F, Edaye S, Huttenhofer A, Brunei C (2007) Anopheles geunbiae miRNAs as 
actors of defence reaction against Plasmodium invasion. NAR 35:6953-6962. 

7. Lcc Y, Kim M, Han J, Yeom KH, Lec S, ct al. (2004) MicroRNA genes arc 
transcribed by RNA polymerase II. The EMBO Journal 23:4051-4060. 

8. Lcc Y, Ahn C, Han J, Choi H, Kim J, ct al. (2003) The nuclear RNasc III 
Drosha initiates microRNA processing. Nature 425:415-419. 

9. Yi R, Qin Y, Macara IG, Cullcn BR (2003) Exportin-5 mediates the nuclear 
export of pre-microRNAs and short hairpin RNAs. Genes and Dev 17:3011- 
3016. 

10. Okada C, Yamashita E, Lcc SJ, Shibata S, Katahira J, et al. (2009) A high- 
resolution structure of the prc-microRNA nuclear export machinery. Science 
326:1275-1279. 

11. Grishok A, Pasquinelli A.E, Contc D, Li N, Parrish S, ct al. (2001) Genes and 
mechanisms related to RNA interference regulate expression of the small 
temporal RNAs that control C.elegans developmental timing. Cell 106:23-3. 

12. Hutvagncr G, McLachlanJ, Pasquinelli AE, Balint E, Tuschl T, et al. (2001) A 
cellular function for the RNA-intcrferencc enzyme Dicer in the maturation of 
the lct-7 small temporal RNA. Science 293:834-838. 

13. Ketting RF, Fischer SEJ, Bernstein E, Sijen T, Hannon GJ, et al. (2001) Dicer 
functions in RNA interference and in synthesis of small RNA involved in 
developmental timing in C. elegans. Genes Dev 15:2654-2659. 

14. Grosshans H, Johnson T, Rcincrt K.L, Gerstein M, Slack FJ (2005) The 
temporal patterning microRNA lct-7 regulates several transcription factors at the 
larval to adult transition in C. elegans. Dev Cell 8:321—330. 

15. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schcltcr JM, ct al. (2005) 
Microarray analysis shows that some MicroRNAs down-regulate large numbers 
of target mRNAs. Nature 433:769-773. 

16. Bahia AC, Dong Y, Blumberg BJ, Mlambo G, Tripathi A, ct al. (2014) 
Exploring Anopheles gut bacteria for Plasmodium blocking activity. Environ 
Microbiol 12381: 1462-2920. 

17. Mathias DK, Pastrana-Mena R, Ranucci E, Tao D, Fcrruti P, et al. (2013) A 
small molecule glyeosaminoglyean mimetic blocks Plasmodium invasion of the 
mosquito midgut. PLoS Pathog. Il:el003757. 

18. Xu P, Vernooy SY, Guo M, Hay BA (2003) The Drosophila MicroRNA Mir- 14 
Suppresses Cell Death and Is Required for Normal Fat Metabolism. Current 
Biol 13:790-795. 

19. Brcnncckc J, Hipfncr DR, Stark A, Russell RB, Cohen SM (2003) bantam 
Encodes a Developmcntally Regulated microRNA that Controls Cell Prolifer- 
ation and Regulates the Proapoptotic Gene hid in Drosophila. Cell 1 13:25—36. 

20. Bejarano F, Smibcrt P, Lai EC (2010) miR-9a prevents apoptosis during wing 
development by repressing Drosophila LIM-only. Dev Biol 338:63-73. 



Table SI Table SI is a table listing miRNAs identified 
containing mismatch with known miRNAs and with An. 
stephensi genome. 

(XLSX) 

Table S2 Table S2 is a table listing miRNAs regulated 
post blood feeding and infection and their targets 
predicted by RNA hybrid tool. 

(XLSX) 

Acknowledgments 

Small RNA Truseq sequencing was performed at NxGenBio Life Sciences, 
New Delhi. We are grateful to the technical assistance of all NIMR staff in 
mosquito collection and parasite infection. Tara Ram's assistance in 
maintenance of mosquitarium at ICGEB is duly acknowledged. 

Author Contributions 

Conceived and designed the experiments: SS RKB. Performed the 
experiments: SJ VR. Analyzed the data: SJ JS. Contributed reagents/ 
materials/ analysis tools: AS SS AT. Wrote the paper: SJ SS RKB. 



21. Lai EC, Tomancak P, Williams RW, Rubin GM (2003) Computational 
identification of Drosophila microRNA genes. Genome Biol 4:42. 

22. Aravin AA, Quintana ML, Yalcin A, Zavolan M, Marks D, ct al. (2003) The 
Small RNA Prolile during Drosophila nulanogaster Development Dev Cell 5:337— 
350. 

23. Ruby JC, Stark A, Johnston WK, Kellis M, Bartel DP, ct al. (2007) Evolution, 
biogenesis, expression, and target predictions of a substantially expanded set of 
Drosophila MicroRNAs. Genome Res 17:1850-1864. 

24. Stark A, Kheradpour P, Parts L, Brcnncckc J, Hodges E, et al. (2007) Systematic 
discovery and characterization of fly microRNAs using 12 Drosophila genomes. 
Genome Research 17:1865-1879. 

25. Mead EA, Tu Z (2008) Cloning, characterization, and expression of microRNAs 
from the Asian malaria mosquito, Anopheles slephensi. BMC Genomics 9:244. 

26. Shrinct J, Jain S, Jain J, Bhatnagar RK, Sunil S (2014) Next Generation 
Sequencing reveals regulation of distinct Aedes microRNAs during chikungunya 
development. PLoS Negl Trop Dis. 8:c2616. 

27. Griffiths-Jones S, Grocock RJ, van Dongcn S, Bateman A, Enright AJ (2006) 
miRBasc: microRNA sequences, targets and gene nomenclature. Nucleic Acids 
Res. 34: 140-144. 

28. Huang Y, Niu B, Gao Y, Fu L, Li W (2010) CD-HIT Suite: a web server for 
clustering and comparing biological sequences. Bioinformatics 26: 680-682. 

29. Kin T, Yamada K, Tcrai G, Okida H, Yoshinari Y, ct al. (2007) IRNAdb: a 
platform for mining/annotating functional RNA candidates from non-coding 
RNA sequences. Nucleic Acids Res 35:D145-148. 

30. Megy K, Emrich SJ, Lawson D, Campbell D, Dialynas E, ct al. (2012) 
VectorBasc: improvements to bioinformatics resource for invertebrate vector 
genomics. Nucleic Acids Res 40: D729-734. 

31. Langmead B, Trapnell C, Pop M, Salzberg SL (2009) Ultrafast and memory- 
efficient alignment of short DNA sequences to the human genome. Genome Biol 
10: R25. 

32. Lorenz R, Bernhart SH, Hcincr Zu Siedcrdissen C, Tafer H, Flamm C, et al. 
(2011) VicnnaRNA Package 2.0. Algorithms Mol Bio 6:26. 

33. Krugcr J, Rchmsmeier M (2006) RNAhybrid: microRNA target prediction easy, 
fast and flexible. Nucleic Acids Res 34: W45 1-454. 

34. Xic C, Mao X, HuangJ, Ding Y, WuJ, ct al. (201 1) KOBAS 2.0: a web server 
for annotation and identification of enriched pathways and diseases. Nucleic 
Acids Res 39: W316-322. 

35. Smoot ME, Ono K, Ruschcinski J, Wang PL, Idckcr T (2011) Cytoscapc 2.8: 
new features for data integration and network visualization. Bioinformatics 27: 
431-432. 

36. Ambros V, Bartel B, Bartel DP, Burge CB, Carrington JC, ct al. (2003) A 
uniform system for microRNA annotation. RNA 9:277-279. 

37. Montalvo Alvarez AM, Landau I, Baccam D (1991) Plasmodium vinckei pcttcri: 
various aspects of its sporogony and exoerythrocytic schizogony. Rev Inst Med 
Trop Sao Paulo 6: 421-6. 

38. Biron DG, Loxdale HD (2013) Host -parasite molecular cross-talk during the 
manipulative process of a host by its parasite. J Exp Biol 216:148-160. 

39. Lefevre T, Thomas F, Schwartz A, Levashina E, Blandin S, ct al. (2007a) 
Malaria Plasmodium agent induces alteration in the head proteome of their 
Anopheles mosquito host. Proteomics 7:1908-1915. 



PLOS ONE | www.plosone.org 



17 



May 2014 | Volume 9 | Issue 5 | e98402 



miRNome Study upon Blood Feeding and Plasmodium Infection in Mosquito 



40. Bryant B, Macdonald W, Raikhcl AS (2010) microRNA miR-275 is 
indispensable for blood digestion and egg development in the mosquito Aedes 
aegypti. PNAS 107:22391-22398. 

41. Hussain M, Walker T, O'Neill SL, Asgari S (2013) Blood meal induced 
microRNA regulates development and immune associated genes in the Dengue 
mosquito vector, Aedes aegypti. Insect Bioehcm Mol Biol 43:146-152. 

42. Bushati N, Stark A, Brcnneckc J, Cohen SM (2008) Temporal reciprocity of 
miRNAs and their targets during the matcrnal-to-zygotic transition in Drosophi/a. 
Curr Biol. 18:50 1-506. 

43. Fullaondo A, Lcc SY (2012) Identification of putative miRNA involved in 
Drosophi/a melanogastcr immune response. Dev Comp Immunol 36:267-273. 

44. Hansen IA, Attardo GM, Park JH, Peng Q, Raikhel AS (2004) Target of 
rapamycin-mediated amino acid signaling in mosquito anautogeny. PNAS 101: 
10626-10631. 

45. Marinotti O, Calvo E, Nguyen QK, Dissanayake S, Ribciro JM, et al. (2006) 
Genome-wide analysis of gene expression in adult Anopheles gambiac. Insect Mol 
Biol. 15:1-12. 

46. Goltsev Y, Rczcndc GL, Vranizan K, Lanzaro G, Valle D, ct al. (2009) 
Developmental and evolutionary basis lor drought tolerance of the Anopheles 
gambiac embryo. Dev Biol. 330:462-70. 

47. Molina-Cruz A, Dcjong RJ, Charles B, Gupta L, Kumar S, ct al. (2007) 
Reactive oxygen species modulate Anopheles gamhiae immunity against bacteria 
and Plasmodium.] Biol Chcm. 6:3217-23. 



48. Goncalvcs RL, Olivcira JH, Oliveira GA, Andersen JF, Olivcira MF, ct al. 
(2012) Mitochondrial reactive oxygen species modulate mosquito susceptibility 
to Plasmodium infection. PLoS One. 7:c41083. 

49. Zugel U, Kaufinann SH (1999) Role of heat shock proteins in protection from 
and pathogenesis of infectious diseases. Clin Microbiol Rev. 12:19-39. 

50. Surachetpong W, Singh N, Cheung KW, Luckhart S (2009) MAPK ERK 
signaling regulates the TGF-bctal -dependent mosquito response to Plasmodium 
falciparum. PLoS Pathog. 5:e 1000366. 

51. Singh CP, Singh J, Nagaraju J (2012) A baculovirus encoded MicroRNA 
(miRNA) suppresses its host miRNA biogenesis by regulating the cxportin-5 
cofactor Ran. J Virol 86: 7867-79. 

52. Li J, Chen Y, Ojn X, Wen J, Ding H, ct al. (2014) miR-138 down-regulates 
miRNA processing in HcLa cells by targeting RMND5A and decreasing 
Exportin-5 stability. Nucleic Acids Res. 42:458-74. 

5.3. Puthiyakunnon S, Yao Y, Li Y, Gu J, Peng H, et al. (2013) Functional 
characterization of three microRNAs of the Asian tiger mosquito, Aedes albopictus. 
Parasit Vectors. 6:2.30. 

54. Triboulct R, Mari B, Lin YL, Chablc-Bessia C, Bcnnasser Y, ct al. (2007) 
Suppression of microRNA-silcncing pathway by HIV-1 during virus replication. 
Science 315:1579-1582. 

55. Navarro L, Jay F, Nomura K, He SY, Voinnct O (2008) Suppression of the 
microRNA pathway by bacterial effector proteins. Science 321:964-967. 



PLOS ONE | www.plosone.org 



18 



May 2014 | Volume 9 | Issue 5 | e98402 



